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ABSTRACT 

Context. The Seyfert 1 galaxy Mrk 841 was observed five times between 2001 and 2005 by the XMM-Newton X-ray observatory. 
The source is well known for showing spectral complexity in the variable iron line and in the soft X-ray excess. 
Aims. The availability of multiple exposures obtained by the reflection grating spectrometer (RGS) cameras allows thorough study of 
the complex absorption and emission spectral features in the soft X-ray band. This paper reports on the first study of Mrk 841 soft 
X-ray spectrum at high spectral resolution. 

Methods. The three combined exposures obtained in January 2001 and the two obtained in January and July 2005 were analysed 
with the SPEX software. 

Results. We detect a two-phase warm absorber. A medium ionisation component (log£~1.5-2.2 ergs s cirr 1 ) is responsible for a deep 
absorption feature in the unresolved transition array of the Fe M-shell and for several absorption lines in the OVI-VIII band, and a 
higher ionisation phase with log£ ~3 ergs s cirr 1 is required to fit absorption in the NelX-X band. The ionisation state and the column 
density of the gas present moderate variation from 2001 to 2005 for both phases. The high ionisation component of the warm absorber 
has no effect on the Fe K band. No significant velocity shift of the absorption lines is measured in the RGS data. Remarkably, the 
2005 spectra show emission features consistent with photoionisation in a high-density (n^lO 11 cirr 3 ) gas. A prominent OVII line 
triplet was clearly observed in January 2005 and narrow radiative recombination continua (RRC) of OVII and CVI were observed in 
both 2005 data sets. A broad Gaussian line around 21.7 A was also required to fit all the data sets. The derived radial distance for the 
emission lines seems to suggest that the photoionisation takes place within the optical broad line region of the source. 

Key words. Galaxies: Seyfert - Galaxies: individuakMrk 841 - Techniques spectroscopic 



1. Introduction 

The model that explains active galactic nuclei (AGN) as being 
powered by accretion on a supermassive black hole surrounded 
by gas (Rees, 1984) is nowadays widely accepted. In the past 
decade, the advent of high-resolution X-ray spectroscopy has al- 
lowed detailed study of ionised material flowing out along the 
line of sight within the central region (e.g. Kaspi et al. 2001; 
Kaastra et al. 2002). Multi-wavelength simultaneous observa- 
tions have shown that X-ray outflows and UV winds are tightly 
connected and that the absorption must take place in a mate- 
rial with multiple ionisation layers (Gabel et al. 2005). More 
than 50% of the Seyfert 1 galaxies are characterised by this gas, 
which is mostly revealed in the X-ray band through detection of 
narrow absorption features from ionised elements, such as car- 
bon, oxygen, neon, iron (Creenshaw et al. 2003, Blustin et al. 
2005). 

The structure of the outflowing wind and its relation to the 
central engine where the AGN continuum radiation originates, 
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are not entirely clear. The absorbing gas seems to be distributed 
diversely in individual sources, its location spanning a scale of 
light days in, e.g., NGC 4051 (Krongold et al. 2007), to a few pc 
in, e.g., NGC 3783 (Behar et al. 2003). 

Variability in the ionised absorbers can often offer a key to 
studying the change in the opacity of the gas in relation to the 
continuum emission. For example, long Chandra observational 
campaigns on bright Seyfert Galaxies have yielded a variety of 
results on the response of the warm absorber to the source lumi- 
nosity fluctuations (Netzer et al. 2002; Krongold et al. 2005). 

More recent studies have shown that the presence of emis- 
sion lines in the soft X-ray band brings additional information 
on the circumnuclear ionised gas. In some Seyfert Is, emission 
lines from highly ionised species can originate very close to or 
within the broad line region of the AGN (Costantini et al. 2007, 
Longinotti et al. 2008), but only in sources with high signal-to- 
noise has it been possible to relate the X-ray emitter with the 
X-ray absorber (NGC 3783, Behar et al. 2003). 

Mrk 841 is a bright Seyfert 1 galaxy at z = 0.0365 (Veron- 
Cetty & Veron, 2001) that has been observed by several X-ray 
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Table 1. XMM-Newton observation log for Mrk 841 



Date 
(yyyy/mm/dd) 


OBSID 


Exp 
(ks) 


1 count rate 
C.ts/s 


(10 cgs) 


2001/01/13 


0112910201 


10 


0.63±0.01 


26.6+1.0 


2001/01/13 


0070740101 


12 


0.72±0.01 


26.6±1.0 


2001/01/14 


0070740301 


14 


0.71+0.01 


26.6±1.0 


2005/01/16 


0205340201 


45 


0.19±0.01 


7.5±0.5 


2005/07/17 


0205340401 


30 


0.23+0.01 


9.1+0.1 



1 In RGS 2 



2 Not corrected for absorption; the 2001 flux refers to the global fit of 3 
data sets 

observatories in the past. Spectral variability was reported in the 
EXOSAT, Ginga and RO SAT data by George et al. (1993), and 
by Nandra et al. (1995). Evidence for a broad iron emission 
line was found in the ASCA spectra reported by Nandra et al. 
(1997). Weaver et al. (2001) have reported marginal evidence for 
line variability. The analysis of broad band BeppoSAX data re- 
vealed a soft X-ray excess and a Compton reflection component 
(Bianchi et al. 2001, 2004), although these authors found that 
warm absorption provided an equally good description of the 
soft X-ray data. Evidence for a moderately deep OVII absorp- 
tion edge was also reported in ASCA data by Reynolds (1997). 
XMM-Newton data at CCD resolution highlighted the presence 
of a variable and complex narrow Fe Ka line (Petrucci et al. 
2002, Longinotti et al. 2004, Petrucci et al. 2007). Besides the 
study of Fe K band, the later paper focused on the analysis of 
the soft excess in Mrk 841 by testing two competing models. 
On one hand, the soft excess could be modelled by reflection 
off a photoionised accretion disc whose spectrum is relativisti- 
cally blurred and smeared (Crummy et al. 2006). A valid alterna- 
tive model was proposed by Gierlinski & Done (2004) who ex- 
plained the soft excess as arising from ionised absorption in a rel- 
ativistically smeared wind. According to Petrucci et al. (2007), 
both models were statistically indistinguishable in the EPIC data 
of Mrk 841. 

This paper presents the very first analysis of the high- 
resolution soft X-ray spectra of Mrk 841, and it makes use of all 
the available XMM-Newton data sets, i.e., the five observations 
previously studied by Petrucci et al. (2007). 

2. The XMM-Newton observations 

Mrk 841 was observed by XMM-Newton 3 times in 2001 and 
twice in 2005 (TableQ]). In this paper, we present only the high- 
resolution spectra gathered by the RGS instrument (den Herder 
et al., 2001). For more details on the EPIC data and on the broad- 
band spectrum, the reader is referred to Petrucci et al. (2007). All 
the raw data were processed with SAS version 7.0.0. Spectral 
files for source and background were produced by running the 
standard procedure rgsproc. Response matrices were created 
by the task rgsrmfgen. The observation performed in January 
2005 presents an anomaly in the illumination of CCD 2 (for both 
RGS1 and RGS2). The background subtraction was checked and 
confirmed as unaffected by this event. 

The source varied significantly in flux and spectral shape 
over the time covered by XMM-Newton observations (see 
Table [TJ and see Petrucci et al. 2007 for an extensive discussion 
on the variability properties of the source). The soft X-ray flux 



decreased by a factor of ~3.5 from 2001 to 2005, thus it is nec- 
essary to consider these two "epochs" separately. 

The three short exposures of 2001 were checked carefully to 
decide whether the six RGS 1 and RGS2 could be combined in 
two spectra (one for each instrument). We fitted the spectra with 
the same model to check the agreement in the spectral shape and 
in the flux level. No significant variation is found in the spectral 
shape in the RGS band, but since the count rates for the 2001 
exposures differ by a few percent (see Table [TJ, instead of co- 
adding the data, we decided to fit simultaneously the six spectra 
of 2001 assuming the same fitting model and normalising the 
data to the RGS 1 in observation 1 1 29 1 020 1 . 

This approach is maintained throughout the paper, and we 
refer to the 2001 simultaneously fitted spectra simply as "Jan 
2001". Since the exposures of 2005 are separated by 6 months, 
we checked whether they could be combined into a single one. 
After fitting the two 2005 spectra with a power law, the spectral 
slopes were found to be different by ~20%, therefore the 2005 
spectra were not co-added in the analysis. 

In summary, the three epoch spectra are kept disjointed dur- 
ing the entire spectral analysis and are referred to as Jan 2001, 
Jan 2005, and July 2005 throughout the paper. 

3. RGS spectral analysis 

The RGS spectral analysis was carried out by using the lat- 
est version of the fitting package SPEX (Kaastra et al. 1996fl 
Galactic absorption of column density N#=2.7x 10 20 cm -2 was 
included in all the following spectral fits (Dickey & Lockman, 
1990). The galactic absorption was modelled with the SPEX HOT 
component with a temperature fixed to 5x1 0~ 4 keV. For all spec- 
tral models, we adopted Solar elemental abundances (Anders & 
Grevesse, 1989). 

The first order spectra were rebinned of a factor of 3 over the 
range 7 to 35 A and the C-statistics was adopted (Cash, 1979). 
Errors are quoted at the 90% confidence level for one interest- 
ing parameter. Throughout the text and tables, the quoted wave- 
lengths have been corrected for the cosmological redshift. 

3.1. The warm absorber 

At first glance, visual inspection of the fluxed spectra (i.e. the 
spectrum plotted in flux units rather than counts/sec) reveals the 
typical signatures of absorption from the so-called unresolved 
transition array of the Fe M shell in the 15-16 A band (e.g. 
Sako et al. 2001). The feature appears to be prominent in Jan 
2005, but it is also easily recognised in the other two data sets. 
We checked the residuals after fitting the spectra with a sim- 
ple power law, confirming absorption in all of them, although 
with different characteristics. This may suggest variability in the 
ionised material along the three epochs and/or the presence of a 
multi-component warm absorber. In order to compare the three 
spectra and test this hypothesis, we took the approach of initially 
fitting them with the same model. 

We constructed a warm absorber model with the SPEX built- 
in model XABS. The ionisation balance was calculated by as- 
suming the spectral energy distribution (SED) of Mrk 841 as the 
ionising continuum spectrum. The SED was constructed with the 
simultaneous optical-UV-to-X-ray fluxes extracted from the op- 
tical monitor and the EPIC instruments onboard XMM-Newton. 
For the 2001 data set the photometric point in the U band (cen- 
tred at 344 nm) was used. In Jan 2005 and July 2005, the source 

1 See also http://www.sron.nl/spex 
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Fig. 1. Spectral energy distribution of Mrk 841 constructed from 
the XMM-Newton data in the optical-X-ray band and from the 
standard AGN radio-IR continuum included in CLOUDY (see 
text for more details). 



NelX (f) OVII RRC2 




Wavelength (A) 

Fig. 2. Jan 2005 RGS 2 data modelled with a single compo- 
nent (dashed red line) and a two-component warm absorber 
(solid blue line). The forbidden emission line of helium-like 
neon and the RRC of OVII, which are included in both models 
(see Sect. 13.21 ). are labelled for clarity. 

was observed with the optical grism, which produces a broader 
band spectrum. The flux was extracted at 344 nm for those data 
sets. At first approximation the complex X-ray unabsorbed broad 
band continuum can be modelled as a broken power law with an 
energy break at about 1 .5 ke V (as derived by a phenomenologi- 
cal fit to the EPIC data, see also Fig. 5 of Petrucci et al. 2007). 
For each of the three data sets, we took the soft power-law in- 
dex from the RGS and the harder energy power law index from 
Petrucci et al. (2007). The very-low energy portion of the SED 
(radio to IR) is the same as the standard AGN continuum used 
in CLOUDY (Ferland et al. 1998). The three SEDs are plotted in 
Fig.rrj Then, we constructed a model for the continuum consist- 
ing of an underlying single power law absorbed by warm ionised 
gas and applied it to the 0.35-1.77 keV data (i.e. the RGS band). 

The column density Nh and the ionisation parameter £ of the 
gas are free parameters. The ionisation parameter is defined as 
L 

£ = -, where L is the ionising luminosity in 1-1000 Rydberg, 

n e R L 

n e the electron density of the gas, and R the distance of the gas 



itself from the central source of radiation. The XABS model ac- 
counts for the width of the absorption lines through the velocity 
dispersion parameter. It has been shown in several Seyfert galax- 
ies that UV and X-ray absorbers originate in the same outflowing 
gas, thus UV observations can in general provide a reliable limit 
on the X-ray line broadening (Arav et al. 2007, Costantini et al. 
2007). Since no constraint on the velocity dispersion of the UV 
lines is readily available in the literature on Mrk 841, this param- 
eter was initially fixed to an average value of 20 km s _1 . When 
it was left free in the fitting process, it tended to decrease to zero 
in all the spectra, thus it was eventually kept fixed to 20 km s . 
Other free parameters in the model are the photon index of the 
power law in the RGS band and the outflowing velocity of the 
warm absorber. The fit statistic improves significantly with the 
addition of the XABS warm absorber model with respect to the 
power-law model, but the fit is still not completely satisfactory 
for all the three epochs. 

A close inspection of the spectra shows that absorption fea- 
tures are visible as negative residuals in the neon band (12- 
14 A). The single warm absorber component is not sufficiently 
ionised to account for absorption in this band. Therefore, a sec- 
ond warm absorber with a higher ionisation level was tested by 
including in the model a second XABS component. The addition 
of the second warm absorber leads to a fit statistic improvement 
by AC=75 and AC=25 for 3 d.o.f. respectively measured in Jan 
and July 2005 spectra. The resulting best-fitting parameters are 
reported in Table [2] The soft photon index refers to the power 
law fitted in the RGS band, i.e. 7-35 A. The quoted parameters 
and C-statistic refer to the model including the emission features 
in Table [3] The improvement is less significant in Jan 2001, but 
this is not surprising since the column densities in Table [2] indi- 
cate that the absorber may be shallower at this epoch. The effect 
of the second warm absorber is shown in Fig.|2]for Jan 2005: the 
Fe UTA is clearly well fitted by the medium ionisation phase, 
whereas the absorption features of NelX ls-3p at 1 1.556 A, and 
NeX Lya at 12.132 A require a more highly ionised plasma to 
be accounted for. 

The two-phase warm absorber provides a good fit to the data 
as can be seen in Fig. [3] The main absorption features are identi- 
fied with transitions of OVI, OVII, OVIII in the 18-23 A band, Fe 
UTA at 15-16 A and Neon IX-X in the 12-14 A. No significant 
velocity shift is detected in any of the two warm absorbers (see 
Table 0, except for the high ionization phase in Jan 2001. The 
upper limits are consistent with outflowing absorbers blueshifted 
by several xlO 2 km s , as commonly observed in Seyfert 1 
Galaxies (e.g. McKernan et al. 2007). 

Even with the addition of the second warm absorber, the fit 
of the Jan 2005 spectrum is not completely satisfactory, since 
positive residuals are still present in the OVII band (21-22 A, see 
top panel of Fig. |4j». The analysis of these residuals in emission 
is reported in the following section. 

3.2. Emission features in the RGS spectra 
3.2.1. Narrow emission lines 

We start by fitting Jan 2005 with three narrow-line components 
at the wavelengths corresponding to the rest frame positions of 
the resonance, intercombination, and forbidden lines (r, i, f) in 
the OVII line triplet at 21.600, 21.790, and 22.101 A, respec- 
tively. The analysis of emission lines in this portion of the spec- 
trum is complex because of the superposition of the absorption 
lines imprinted by the warm absorber. The peak wavelengths of 
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Table 2. Best-fitting parameters for the two-phase warm absorber. 



Obs 


r« 


ft 


Log£ 
(ergs s cm 1 ) 


N H 
10 21 cm" 2 


V 

km s -1 


Log^ 
(ergs s cm 1 ) 


10 21 cm" 2 


V 

km s" 1 


Cstat/d.o.f. 


Jan 01 
Jan 05 
July 05 


2.60+ 
2.4U 
2.07i 


0.07 
.0.07 
-0.07 


2.2±0.2 
1.5+0.1 

1 7 +a3 

1 -'-0.2 


1.2+0.5 
2.3±0.5 

3 9+ 2 - 6 


<1600 
< 100 
<1000 


T Q+0.8 

J -0.4 

3.2±0.1 
2.8±0.3 


<17.0 
30 +uo 

■'"■"-20.0 

7 6 +ls 

'• u -5.5 


1000+460 
<100 
<1000 


3082/2725 
1958/1435 
1855/1437 




Jan 2001 




0.0000 
0.0006 

V to 0.0004 

CM 
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o 

§ 0.0002 

o 
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0.0000 
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Fig. 3. From top to bottom: RGS rest-frame spectra of Jan 2001 (for plotting purpose only obsid 0070740101 is shown), Jan 2005 
and July 2005 fitted with the best-fitting model including the components in Tables [2] and [3] For plotting purposes, the spectra are 
re-binned by a factor of 9 and plotted in the 7-30 A band, where most of the features of interest fall. The labels mark the ionic 
species with the highest contribution in our warm absorber model. 



the three lines are kept frozen to the laboratory values. Since 
the r transition has a large oscillator strength in the absorption 
component, most of the emission line at this wavelength will be 
re-absorbed, making a precise measure of its flux somewhat dif- 
ficult. The i and r components are detected, while only an upper 
limit is found for the/line. Less prominent positive residuals are 
present in the NelX band, hence we searched for emission lines 
also in this portion of the spectrum. Only the NelX forbidden 
component at 13.699 A is detected (the neon band is shown in 
Fig. The emission lines fluxes and significances in terms of 
C statistic improvement for one free parameter, are reported in 
Table[3] The structure of the line triplet fitted in Jan 2005 data is 
plotted in Figure [4] (thick red line). 

We model the spectra of Jan 2001 and July 2005 with the 
same emission line components included for Jan 2005. Not sur- 



prisingly, given that the continuum flux is higher in these two 
data sets (see Table [TJ, we could measure only the upper limits 
on the line fluxes, which are generally consistent with Jan 2005 
spectra, except for the OVII resonance line in July 2005 that may 
be severely underestimated. 



3.2.2. Radiative recombination continua 

The spectrum of Jan 2005 displays an additional emission fea- 
ture around 17 A, consistent with the position of the radiative 
recombination continuum (RRC) of OVII (see Fig.|2]for a close- 
up of this feature). A closer look at longer wavelengths suggests 
that other RRC components from Carbon may give significant 
contribution in the positive residuals of the data. In our analysis, 
these features are fitted with the RRC model in SPEX. The width 
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3.2.3. Broad oxygen line 

Once the narrow lines are included in the fit, the residu- 
als around the OVII triplet in Jan 2005 still show some 
structure in excess, particularly on the left side of the reso- 
nant line. To account for these residuals, a broad Gaussian 
line was added in the model with width and wavelength 
free to vary. The resulting best-fitting parameters for this 
component are A=21.75+°'^ A, FWHM=0.60+°g A, and 

flux=0.64+0.53xl0~ 4 photons cirT 2 s -1 . The addition of the 
broad line component leads to an improvement of AC stat= 1 8 
for 3 degrees of freedom, corresponding to a significance level 
higher than 99%. Fig. [4] shows the separate contribution of the 
Gaussian profile to the best fit model for the OVII line triplet. 
For simplicity, the broad line at 2 1 .7 A will only be referred to as 
"OVII line", without identifying a specific transition. The same 
line component was included in July 2005 yielding the following 



best-fitting parameters: ,1=21.71 ±0.23 A, FWHM=0.54+°jI A 
and flux=0.51+0.38xl0- 4 photons cirT 2 



20.5 



21.0 



21.5 22.0 
Wavelength (A) 



22.5 



23.0 



0.44 

s ' . Although the sig- 
nificance of this component in July 2005 (around 90%) is lower 
than in Jan 2005, the generally good agreement in the best- 
fitting parameters for the two spectra seem to support the pres- 
ence of the same line in July 2005. In Jan 2001 data, we added 
the broad component as we did for the narrow lines. Since the 
higher continuum flux prevents any determination of the line 
properties, to measure the line flux, the position and width of the 
Gaussian component have been fixed to the average value found 
for the 2005 epochs (see Table[3]l. An upper limit of <0.85 pho- 
tons irT 2 s _1 (consistent with the 2005 measurements) was found 
for the broad component in 2001 . 



Fig. 4. Top: Residuals of the continuum model in the OVII band. 
Bottom: Zoom on the OVII triplet in Jan 2005. Line parame- 
ters are reported in Table 3. The solid red line corresponds to 
the best-fitting model (including the emission lines in Table [3] 
and the underlying warm absorber in Tabled. The dashed blue 
line shows the underlying OVII broad Gaussian line after remov- 
ing the contribution of the narrow lines triplet and the warm ab- 
sorber. 



of the RRC profiles is a very sensitive indicator of the electron 
temperature in the recombining plasma (Liedahl & Paerels 1996, 
Liedahl 1999). Our data suggest the presence of a narrow RRC 
at the OVII wavelength (see Fig. |5J, indicative of temperature 
of fewxlO 4 K. Therefore, we fit the RRC emission measur^ 
of OVII and CVI as free parameters assuming an electron tem- 
perature corresponding to kT- 3 eV, i.e. T=3.5xl0 4 K. If the 
temperature is left free to vary in the fit, we find kT=2.8^^ eV. 

The improvement in the fit statistic for adding the RRC com- 
ponent in Jan 2005 is AC=23 for 2 degrees of freedom. The 
OVII RRC seems to be slightly redshifted from its laboratory 
position; we applied a small redshift of z=0.001 to match the ob- 
served excess in the data. The same was done for the C VI RRC 
around 25 A, where the shift is less evident due to the lower 
signal-to-noise of our spectra in this band. As described in the 
previous section for the narrow emission lines, we included the 
RRC in the model for all the three spectra and reported the best- 
fitting parameters in Table [3] 



2 The emission measure corresponds to n e nf on V, where n e is the gas 
electron density, n, on the density of the parent ion that recombines to the 
ground state, and V the emitting volume 



4. Discussion 

The analysis of the RGS data of Mrk 841 presented herein pro- 
vides a considerable amount of information. A large contribution 
from a warm absorber consisting of two ionisation components 
has been revealed and most interestingly, it proves to be rela- 
tively stable on a time scale of 4 years. In addition, emission 
features from highly ionised species of neon, oxygen, and car- 
bon have been ascertained. These topics will be addressed sep- 
arately in the following sections. The discussion is based on the 
results from the global best fit (i.e. absorption+emission), which 
is reported in Tables [2] and [3] 

4.1. A two-phase warm absorber in Mrk 841 

The properties of the ionised absorber emerging from the spec- 
tral fits described in Sect. [XT] indicate that there are two distinct 
levels of ionisation in the absorbing gas. A medium-phase warm 
absorber characterised by log£ in the range 1.5-2.2 ergs s crrT 1 
and column density of the order of 10 21 crrT 2 is responsible for 
transition in the Fe UTA and in the oxygen band. 

The availability of spectra from different epochs may allow 
a variability study of the absorber physical conditions, particu- 
larly considering that because the Fe UTA is the result of a blend 
of absorption lines arising from several iron ions, its shape and 
position are very sensitive to changes in the ionisation balance 
of the gas and, ultimately, to luminosity variations of the source 
(Behar et al., 2001). The errors in Table|2]show that £ and N# are 
generally consistent through the three data sets, although there 
seems to be a decrease in ionisation accompanied by a moder- 
ate increase of the column density going from 2001 to Jan 2005. 
Indeed, the Fe UTA is much more evident in both 2005 spectra 
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Table 3. Emission lines parameters and emission measure for RRC in the three epochs. 



Obs 


Line 


,1 

(A) 


Flux 

(10- 4 phcnr 2 s" 1 ) 


AC 1 


FWHM 

(A) 


RRC EM 

(10 5S cm" 3 ) 


Jan 01 


OVn (r) 
OVII (i) 
OVII (f) 
NelX (f) 


21.600 
21.790 
22.101 
13.699 


<0.42 
< 0.15 
<0.43 
<0.12 


- 


0-width 
0-width 
0-width 
0-width 


OVII <450 
CVR4300 




OVU (gau) 


21.7 


<0.85 


- 


0.55 




Jan 05 


OVII (r) 
OVII (i) 
OVII (f) 
NelX (f) 


21.600 
21.790 
22.101 
13.699 


0.24±0.16 
26 +0 - 26 
<0.20 
0.10±0.05 


17 
16 

10 


- 


OVII= 1500+800 
CVI=2500±1800 
AC stat=23 




OVII (gau) 


21 75 +0 - 23 


0.64±0.53 


18 


60 +0 52 

"•""-0.35 




July 05 


OVII (r) 
OVII (i) 
OVII (f) 
NelX (f) 


21.600 
21.790 
22.101 
13.699 


<0.04 
<0.25 
<0.23 
<0.15 






OVII=1860±1000 
CVI=2500±2300 
AC stat=15 




OVII (gau) 


21.71±0.23 


0.51±0.38 


6 


54 +a57 





'When AC stat <3, the upper limits on the line flux are reported 



because the lower ionisation state observed at this epoch is able 
to produce a high number of transitions in Fe XV-XVII in the 
15-16 A region. According to our model, absorption from OVI 
and OVII is also strong in this source. These absorption lines are 
not easily recognisable, especially in Jan 2005 spectrum, due to 
the presence of the OVII emission line triplet. The OVI-OVII ab- 
sorption is likely associated to the gas component in a medium 
ionisation state. 



The higher ionisation phase (log^~3 ergs s cm l ) is instead 
responsible for imprinting absorption at shorter wavelengths 
mainly in the neon band around 12-14 A. This component is 
characterised by a column density of the order of 10 22 cm -2 , 
much higher than the one measured in the medium phase. In 
principle, warm gas at this ionisation level can produce absorp- 
tion lines in the Fe K band around 6-7 keV, which may play 
a key role in shaping the Fe K alpha line profile (e.g. Revees 
et al. 2004, Turner & Miller, 2009). Indeed, Longinotti et al. 
(2004) and Petrucci et al. (2007) reported peculiar behaviour in 
the Fe Ka line in their previous analysis of Mrk 841, but none of 
these works included detailed modelling of the warm absorber. 
The detection of a high ionisation warm absorber in the RGS 
data prompted us to test for absorption effects on the Fe K line. 
This check was performed by simply applying our best-fitting 
model to the EPIC-pn data of Jan 2005, i.e., the spectrum with 
the highest column density in the high ionisation phase (Table|2]). 
With the parameters frozen to the values of Table |2l we found 
that this warm absorber component is not deep enough for intro- 
ducing any additional curvature around the Fe line, or for modi- 
fying the line profile itself. 



4.2. Density of the warm absorber 

In about four-years time, the low-energy flux of Mrk 841 de- 
creased by a factor of 3.5 (Table 1). In Jan 2005, when the source 
was fainter, we measured also a lower ionisation parameter for 
the cold phase, while the column density stayed constant within 
the errors (Table 2). Given the uncertainty on the outflow veloc- 
ity, we assume here that we detected the same gas material in 
both 2001 and 2005, albeit with different ionisations. 

We estimated a lower limit for the density of such gas, as- 
suming a recombination time of 1422 days (4 years time between 

to 

Jan 2001 and Jan 2005), i.e. t s =- in the source rest frame. 



1 +z 



We considered all the Fe ions detected in the spectrum. Each ion 
has its own recombination time, which is inversely related to the 
gas density n e , as reported in Bottorff et al. (2000) and Detmers 
et al. (2008): 

f(X M ) ariXt-df- 1 



t{X i ) = \a r {X i )n e 



f(Xi) a r (Xd 
where f(X,) is the recombination timescale of the ion Xj, 
or(Xj) is the recombination rate from ion Xj = \ to ion Xj, and f(Xj) 
represents the fraction of the element X at the ionisation state i. 

This relation allows us to put further constraints on the lower 
limit of the density. We used the values predicted by the best-fit 
global modelling of the absorber, as the quality of the data does 
not allow a line-by-line fitting on the Fe ions to be performed, 
especially in the crowded region of the iron UTA. The associated 
error on the ionic column densities is therefore the scaled value 
of the error on the total Nh- 

The gas density is a function of the recombination rates, 
which are in turn a function of the gas temperature T e , before 
recombination (Bottorff et al. 2000, Detmers et al. 2008). The 
value of T e is linked to the ionisation parameter £ via the ion- 
isation balance. The error on £ (Table 2), rather than the error 
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Fig. 5. Electron density in the warm absorber estimated using 
the recombination of the iron ions produced by the gas with 
log£ = 2.2 ± 0.23 and 7Y H ~ 1.2 ± 0.5 x 10 21 cirr 2 . 



on individual ionic column densities, has the major weight in 
the density determination. The radiative recombination rates for 
iron were analytically determined using the expression and re- 
combination coefficient provided in Woods et al. (1981). 

In Fig. [5] we show the density lower limits for the Fe ions. 
The ion that provides the higher value for n = 3. 1 ± 1 .4x 10 3 cirr 3 
is Fe xvn. 



4.3. On the origin of the narrow emission lines and radiative 
recombination continua 

The following discussion is based chiefly on the findings 
for Jan 2005 spectrum, where the emission line components 
are detected with highest significance. The narrowness of the 
RRC profiles is consistent with an electron temperature of a 
few xlO 4 K, and it seems to indicate that photoionisation from 
the central source is the primary mechanism acting on the sur- 
rounding gas (e.g. Kinkhabwala et al. 2002, Armentrout et al. 
2007). 

The relative intensity of line triplets in He-like ions provides 
an excellent diagnostic of the physical properties of the emit- 
ting plasma (Porquet & Dubau, 2000 and reference therein). In 
general, it has been observed that photoionisation processes in 
Seyfert Galaxies give rise to line triplets dominated by the for- 
bidden transition (e.g. NGC 1068, Kinkhabwala et al. 2002). If 
we consider the line fluxes of the OVII triplet in Mrk 841, we 
note that their ratio is not typical of photoionisation, the less 
prominent line being the forbidden transition, as is also evident 
from Fig. [4] The intercombination component at 21.790 A is as 
strong as the resonance one, if not higher: the presence of a bad 
pixel between the r and i lines does not allow precise reconstruc- 
tion of the r emission line. On the contrary, we measure an upper 
limit for the forbidden component. Since this line is not affected 
by instrumental effects and is generally expected to be dominant 
in AGN spectra, this measurement provides a stringent limit to 

be considered in the discussion of the line ratios. 

/ 

The R ratio - provides a compelling probe of the gas den- 
i 

sity in photoionised plasma (see Porquet & Dubau, 2000). The 
low value of this ratio (R <0.76) in Mrk 841 points to emission 
from a gas at very high density, i.e. n e > 10 11 cirr 3 . This result 



partly recalls the case of the Seyfert 1 Galaxy Mrk 335 that dis- 
played a photoionised OVII triplet with a dominant i line during 
a particularly low flux state (Longinotti et al. 2008). 



The G ratio 



f + i 



instead is often used as an estimator of 



the electron temperature in the gas. We can estimate an upper 
limit on the G ratio by taking the upper limits on the three com- 
ponents of the triplet into account (i.e., the upper error bars for 
the two lines for which they are available and the upper limit on 
the /line). In this way, the G ratio is estimated to be less than 
1.8. This limit is somewhat intermediate between a pure colli - 
sional and a pure photoionisation-dominated plasma (see Fig. 7 
in Porquet & Dubau, 2000), possibly indicating that mixed pro- 
cesses are at work. We note some discrepancy between the tem- 
perature of the recombining plasma inferred in this work by the 
RRC width (~10 4 K) and the temperature prescribed for the hy- 
brid plasma by the calculations reported in Porquet & Dubau 
(2000). However, as highlighted for Mrk 3 (Sako et al. 2000) and 
NGC 4151 (Armentrout et al. 2007), the intensity of resonance 
lines in He-like triplets could be enhanced by photo-excitation 
mechanisms induced by the AGN, which leads to a decrease in 
the observed G ratio (Godet et al. 2004). In principle, this effect 
could be tested in our data by accounting for microturbulence 
in the resonance line (see Armentrout et al. 2007), but we recall 
that the velocity dispersion was undetermined even for the warm 
absorber lines (less than 20kms~ 1 , see Sect. 13. lb . thus not easily 
determined with the present instrumental resolution. 

4.4. One photoionisation region for the warm absorber and 
the warm mirror? 

While certainly very interesting, the soft X-ray properties in 
Mrk 841 are not unique because several AGNs show fea- 
tures from ionised plasma both in absorption and emission e.g. 
NGC 4151, (Armentrout et al. 2007, Kraemer et al. 2005) and 
NGC 3783 (Behar et al. 2003). For some cases, attempts are be- 
ing made to ascribe the two processes to the same photoionised 
plasma (e.g. NGC 4051, Nucita et al. submitted). The main un- 
known is the geometry of the circumnuclear gas, which should 
act as a "warm mirror" for scattering the emission lines and, at 
the same time, it should intercept the observer's line of sight for 
producing warm absorber features. 

For a large number of X-ray obscured Seyfert galaxies, there 
are several indications that the X-ray photoionisation region 
where the bulk of soft X-ray line emission takes place, is spa- 
tially coincident with the optical narrow line region (NLR) (see 
Bianchi et al. 2006 and Guainazzi & Bianchi 2007). Seyfert 1 
objects, with unobscured X-ray spectra, offer instead a direct 
view of the continuum source and the surrounding high density 
gas giving rise to optical broad lines (i.e. the broad line region 
(BLR)), therefore introducing more diversity in determining the 
distribution of the hot plasma. As an example, in the Seyfert 1 
NGC 4051 (Pounds et al. 2004) and NGC 5548 (Kaastra et al. 
2002, Detmers et al. 2008), the gas density inferred from the ob- 
served emission lines indicates an origin in the NLR, whereas 
in Mrk 335, the high density line-emitting gas (n<,=10 9 ~ n cm -3 ) 
was located within the source BLR (Longinotti et al. 2008). 

The present study of Mrk 841 with the detection of an im- 
portant warm absorber component and photoionised emission 
lines, provides the opportunity to compare and, possibly, link 
their properties. The ionisation parameter of the photoionised 
L 



gas 



n„R 2 



depends on the plasma electron density, its radial 



distance from the source of radiation, and the ionising lumi- 
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nosity. An average luminosity of 10 ergs s was estimated 
from the SED in Fig. [T]by integrating the data of Jan 2001 and 
Jan 2005 over 1-1000 Rydberg. 

In Sect. 14.21 we derived a lower limit for the density of the 
warm absorber (> 10 3 cm -3 ), which can loosely constrain the dis- 
tance to < few tens of pc. Nevertheless, this constraint can be 
refined by using the limit on the density inferred by the OVII 
triplet (n e > 10 11 cm -3 ), under the assumption that the emission 
lines originate within a gas at the same ionisation state of the 
warm absorber gas. Our two-phase model for Jan 2005 yields 
log£~1.5 and 3.2 ergs s cm -1 (see Table 0. However, in the 
model of the gas with log£~3.2 ergs s cm -1 the column density 
of OIX is much higher than the one of OVII and OVIII, therefore 
such a highly ionised gas is unlikely to produce OVII emission. 
Assuming the ionisation parameter of the low phase, the radial 
distance R of the "warm emitter gas" is then constrained to less 
than 1.6xl0 16 cm. 

The person who wishes to put this result in the context of 
the source may consider that the optical BLR lies at a radius 
of about 4.4xl0 16 cm, as roughly estimated from the FWHM 
of 5500 km s -1 in the optical lines (Boroson & Green, 1992). 
Thus, the ionised gas seems to encompass the BLR, analogously 
to Mrk 335 (Longinotti et al. 2008). We checked the agreement 
of the optical FWHM and the width of the OVII line triplet by 
fitting the three components with a Gaussian profile and freez- 
ing the rest of the model. We found the following FWHM for 
r, i, f, respectively: 400+60, <1400 and <5300 km s _1 . While 
the velocities in the i and / components can still be reconciled 
with BLR, it is hard to explain the apparent inconsistency of the 
resonant line with the rest of the triplet. We note, however, that 
the reconstruction of the r line is rather complex, as highlighted 
in Sect. 13.2.11 A more sophisticated modelling of the emission 
lines in high signal-to-noise data would be needed to supplement 
our hypothesis. 

A few words must be said on the line fluxes in the other two 
data sets for which only upper limits are available. The appar- 
ent lack of strong emission lines in Jan 2001 and Jan 2005 can 
be phenomenologically explained as an effect of the increased 
continuum flux. As a consequence, the emission component is 
swamped and covered in the observed spectra. While it is hard 
to exclude intrinsic variability in the emission line fluxes, espe- 
cially since the continuum flux varies, there is no compelling 
evidence for variation in the present data. 

4.5. The broad OVII line in the RGS data 

The broad emission line around ~21.7 A is detected with high 
significance in both 2005 observations. If fitted with a Gaussian 
profile, the FWHM corresponds to ~8200!™» km s _1 , for 
Jan 2005. Considering the large error bars, this velocity is con- 
sistent with the optical H/3 FWHM of 5500 km s -1 reported by 
Boroson & Green (1992). Broad emission features in the soft 
X-ray spectra of Seyfert 1 Galaxies has been reported in the 
Chandra-LETG of Mrk 279 (Costantini et al. 2007), NGC 5548 
(Steenbrugge et al. 2005), NGC 4051 (Steenbrugge et al. 2009), 
and in the XMM-Newton-RGS spectra of NGC 405 1 (Ogle et al. 
2004, Ponti et al. 2006). In most of these cases, the line profiles 
were consistent with Gaussians, and in the first two sources, the 
velocity width of the emission lines was in agreement with mea- 
surements of UV lines in the BLR. 

In the past, it was proposed to model RGS spectral lines 
from H-like ions in the Seyfert 1 Galaxies of MCG-6-30-15 and 
Mrk 766 with a relativistic profile (Branduardi-Raymont et al. 
2001). We tried to convolve the OVII Gaussian emission line in 



Mrk 841 with a relativistic profile using the LAOR kernel (Laor, 
1991) in Jan 2005 spectrum. No good fit was found, perhaps 
not surprisingly since the data quality of Mrk 841 is not com- 
parable to the long-integration RGS spectra of the two bright 
sources mentioned above, keeping the test for relativistic broad- 
ening from being viable in this source. 

Considering the agreement of the OVII and H/3 lines FWHM 
in Mrk 841, a BLR origin for the OVII broad line is favoured. 
Unfortunately, the lack of simultaneous UV spectral data for this 
source prevents us from performing a more detailed analysis of 
this feature and its relation to the multi-wavelength properties. 

5. Summary of conclusions 

In the following we summarise the main results of the first anal- 
ysis of Mrk 841 high-resolution X-ray spectra: 

- Warm absorber: a two-phase warm absorber with log£~1.5- 
2.2 ergs s cm -1 and log£~3 ergs s cm -1 was observed in the 
three spectra from 2001 through 2005, a moderate decrease 
in the ionisation state was observed going from the brightest 
flux state (2001) to the dimmer state (2005), no outflow was 
detected in the data and the high ionisation component has 
no effect on the Fe K band. 

- Emission features: narrow OVII and NelX emission lines 
were detected in the Jan 2005 spectrum, the line ratio in the 
OVII triplet points to an origin in a gas at high electron den- 
sity (n c >10 n cm -3 ), the presence of narrow RRC from OVII 
and CVI is consistent with photoionisation processes, if it 
is assumed that the emission component and the warm ab- 
sorber originate within a gas at the same ionisation state, the 
inferred distance for the emission lines is constrained within 
10 16 cm from the central source. 

- Broad Gaussian: a broad Gaussian line centred at 21.7 A 
was interpreted as emission from OVII, the FWHM is fully 
compatible with that of the optical broad emission lines in 
Mrk 841. 
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